Genome-scale Metabolic Network Research Articles

Spatial transcriptome-guided multi-scale framework connects P. aeruginosa metabolic states to oxidative stress biofilm microenvironment.

With the generation of spatially resolved transcriptomics of microbial biofilms, computational tools can be used to integrate this data to elucidate the multi-scale mechanisms controlling heterogeneous biofilm metabolism. This work presents a Multi-scale model of Metabolism In Cellular Systems (MiMICS) which is a computational framework that couples a genome-scale metabolic network reconstruction (GENRE) with Hybrid Automata Library (HAL), an existing agent-based model and reaction-diffusion model platform. A key feature of MiMICS is the ability to incorporate multiple -omics-guided metabolic models, which can represent unique metabolic states that yield different metabolic parameter values passed to the extracellular models. We used MiMICS to simulate Pseudomonas aeruginosa regulation of denitrification and oxidative stress metabolism in hypoxic and nitric oxide (NO) biofilm microenvironments. Integration of P. aeruginosa PA14 biofilm spatial transcriptomic data into a P. aeruginosa PA14 GENRE generated four PA14 metabolic model states that were input into MiMICS. Characteristic of aerobic, denitrification, and oxidative stress metabolism, the four metabolic model states predicted different oxygen, nitrate, and NO exchange fluxes that were passed as inputs to update the agent's local metabolite concentrations in the extracellular reaction-diffusion model. Individual bacterial agents chose a PA14 metabolic model state based on a combination of stochastic rules, and agents sensing local oxygen and NO. Transcriptome-guided MiMICS predictions suggested microscale denitrification and oxidative stress metabolic heterogeneity emerged due to local variability in the NO biofilm microenvironment. MiMICS accurately predicted the biofilm's spatial relationships between denitrification, oxidative stress, and central carbon metabolism. As simulated cells responded to extracellular NO, MiMICS revealed dynamics of cell populations heterogeneously upregulating reactions in the denitrification pathway, which may function to maintain NO levels within non-toxic ranges. We demonstrated that MiMICS is a valuable computational tool to incorporate multiple -omics-guided metabolic models to mechanistically map heterogeneous microbial metabolic states to the biofilm microenvironment.

Metabolic modeling of Halomonas campaniensis improves polyhydroxybutyrate production under nitrogen limitation

Poly-hydroxybutyrate (PHB) is an environmentally friendly alternative for conventional fossil fuel-based plastics that is produced by various microorganisms. Large-scale PHB production is challenging due to the comparatively higher biomanufacturing costs. A PHB overproducer is the haloalkaliphilic bacterium Halomonas campaniensis, which has low nutritional requirements and can grow in cultures with high salt concentrations, rendering it resistant to contamination. Despite its virtues, the metabolic capabilities of H. campaniensis as well as the limitations hindering higher PHB production remain poorly studied. To address this limitation, we present HaloGEM, the first high-quality genome-scale metabolic network reconstruction, which encompasses 888 genes, 1528 reactions (1257 gene-associated), and 1274 metabolites. HaloGEM not only displays excellent agreement with previous growth data and experiments from this study, but it also revealed nitrogen as a limiting nutrient when growing aerobically under high salt concentrations using glucose as carbon source. Among different nitrogen source mixtures for optimal growth, HaloGEM predicted glutamate and arginine as a promising mixture producing increases of 54.2% and 153.4% in the biomass yield and PHB titer, respectively. Furthermore, the model was used to predict genetic interventions for increasing PHB yield, which were consistent with the rationale of previously reported strategies. Overall, the presented reconstruction advances our understanding of the metabolic capabilities of H. campaniensis for rationally engineering this next-generation industrial biotechnology platform.Key pointsA comprehensive genome-scale metabolic reconstruction of H. campaniensis was developed.Experiments and simulations predict N limitation in minimal media under aerobiosis.In silico media design increased experimental biomass yield and PHB titer.

Reconstruction of Metabolic-Protein Interaction Integrated Network of Eriocheir sinensis and Analysis of Ecdysone Synthesis.

Integrated networks have become a new interest in genome-scale network research due to their ability to comprehensively reflect and analyze the molecular processes in cells. Currently, none of the integrated networks have been reported for higher organisms. Eriocheir sinensis is a typical aquatic animal that grows through ecdysis. Ecdysone has been identified to be a crucial regulator of ecdysis, but the influence factors and regulatory mechanisms of ecdysone synthesis in E. sinensis are still unclear. In this work, the genome-scale metabolic network and protein-protein interaction network of E. sinensis were integrated to reconstruct a metabolic-protein interaction integrated network (MPIN). The MPIN was used to analyze the influence factors of ecdysone synthesis through flux variation analysis. In total, 236 integrated reactions (IRs) were found to influence the ecdysone synthesis of which 16 IRs had a significant impact. These IRs constitute three ecdysone synthesis routes. It is found that there might be alternative pathways to obtain cholesterol for ecdysone synthesis in E. sinensis instead of absorbing it directly from the feeds. The MPIN reconstructed in this work is the first integrated network for higher organisms. The analysis based on the MPIN supplies important information for the mechanism analysis of ecdysone synthesis in E. sinensis.

Reconstruction of the genome-scale metabolic network model of Sinorhizobium fredii CCBAU45436 for free-living and symbiotic states.

Sinorhizobium fredii CCBAU45436 is an excellent rhizobium that plays an important role in agricultural production. However, there still needs more comprehensive understanding of the metabolic system of S. fredii CCBAU45436, which hinders its application in agriculture. Therefore, based on the first-generation metabolic model iCC541 we developed a new genome-scale metabolic model iAQY970, which contains 970 genes, 1,052 reactions, 942 metabolites and is scored 89% in the MEMOTE test. Cell growth phenotype predicted by iAQY970 is 81.7% consistent with the experimental data. The results of mapping the proteome data under free-living and symbiosis conditions to the model showed that the biomass production rate in the logarithmic phase was faster than that in the stable phase, and the nitrogen fixation efficiency of rhizobia parasitized in cultivated soybean was higher than that in wild-type soybean, which was consistent with the actual situation. In the symbiotic condition, there are 184 genes that would affect growth, of which 94 are essential; In the free-living condition, there are 143 genes that influence growth, of which 78 are essential. Among them, 86 of the 94 essential genes in the symbiotic condition were consistent with the prediction of iCC541, and 44 essential genes were confirmed by literature information; meanwhile, 30 genes were identified by DEG and 33 genes were identified by Geptop. In addition, we extracted four key nitrogen fixation modules from the model and predicted that sulfite reductase (EC 1.8.7.1) and nitrogenase (EC 1.18.6.1) as the target enzymes to enhance nitrogen fixation by MOMA, which provided a potential focus for strain optimization. Through the comprehensive metabolic model, we can better understand the metabolic capabilities of S. fredii CCBAU45436 and make full use of it in the future.

A dynamic model of growth phase of bio-conversion of methane to polyhydroxybutyrate using dynamic flux balance analysis.

Biological conversion of waste methane to biodegradable plastics is a way of reducing their production cost. This study addresses the computational modeling of the growth phase reactor of the process of polyhydroxybutyrate production. The model was used for investigating the effect of gas recycling and inlet gas retention time on the reactor performance. The model was run by the use of a genome-scale metabolic network of Methylocystis hirsuta in a dynamic flux balance analysis framework. The reactor has been modeled for two separate feeding scenarios: a pure methane feed and a biogas feed. The mass transfer coefficient parameter was predicted as a function of superficial gas velocities by the regression of data from published experiments. The results show an increase of removal efficiency by 38% and biomass concentration by 2.8g/L with the increase of gas recycle ratio from 0 to 30 at the empty bed residence time of 60 .

Identifying metabolic adaptations characteristic of cardiotoxicity using paired transcriptomics and metabolomics data integrated with a computational model of heart metabolism.

Improvements in the diagnosis and treatment of cancer have revealed long-term side effects of chemotherapeutics, particularly cardiotoxicity. Here, we present paired transcriptomics and metabolomics data characterizing in vitro cardiotoxicity to three compounds: 5-fluorouracil, acetaminophen, and doxorubicin. Standard gene enrichment and metabolomics approaches identify some commonly affected pathways and metabolites but are not able to readily identify metabolic adaptations in response to cardiotoxicity. The paired data was integrated with a genome-scale metabolic network reconstruction of the heart to identify shifted metabolic functions, unique metabolic reactions, and changes in flux in metabolic reactions in response to these compounds. Using this approach, we confirm previously seen changes in the p53 pathway by doxorubicin and RNA synthesis by 5-fluorouracil, we find evidence for an increase in phospholipid metabolism in response to acetaminophen, and we see a shift in central carbon metabolism suggesting an increase in metabolic demand after treatment with doxorubicin and 5-fluorouracil.

Genome scale metabolic network modelling for metabolic profile predictions.

Metabolic profiling (metabolomics) aims at measuring small molecules (metabolites) in complex samples like blood or urine for human health studies. While biomarker-based assessment often relies on a single molecule, metabolic profiling combines several metabolites to create a more complex and more specific fingerprint of the disease. However, in contrast to genomics, there is no unique metabolomics setup able to measure the entire metabolome. This challenge leads to tedious and resource consuming preliminary studies to be able to design the right metabolomics experiment. In that context, computer assisted metabolic profiling can be of strong added value to design metabolomics studies more quickly and efficiently. We propose a constraint-based modelling approach which predicts in silico profiles of metabolites that are more likely to be differentially abundant under a given metabolic perturbation (e.g. due to a genetic disease), using flux simulation. In genome-scale metabolic networks, the fluxes of exchange reactions, also known as the flow of metabolites through their external transport reactions, can be simulated and compared between control and disease conditions in order to calculate changes in metabolite import and export. These import/export flux differences would be expected to induce changes in circulating biofluid levels of those metabolites, which can then be interpreted as potential biomarkers or metabolites of interest. In this study, we present SAMBA (SAMpling Biomarker Analysis), an approach which simulates fluxes in exchange reactions following a metabolic perturbation using random sampling, compares the simulated flux distributions between the baseline and modulated conditions, and ranks predicted differentially exchanged metabolites as potential biomarkers for the perturbation. We show that there is a good fit between simulated metabolic exchange profiles and experimental differential metabolites detected in plasma, such as patient data from the disease database OMIM, and metabolic trait-SNP associations found in mGWAS studies. These biomarker recommendations can provide insight into the underlying mechanism or metabolic pathway perturbation lying behind observed metabolite differential abundances, and suggest new metabolites as potential avenues for further experimental analyses.

Machine learning predicts system-wide metabolic flux control in cyanobacteria

Metabolic fluxes and their control mechanisms are fundamental in cellular metabolism, offering insights for the study of biological systems and biotechnological applications. However, quantitative and predictive understanding of controlling biochemical reactions in microbial cell factories, especially at the system level, is limited. In this work, we present ARCTICA, a computational framework that integrates constraint-based modelling with machine learning tools to address this challenge. Using the model cyanobacterium Synechocystis sp. PCC 6803 as chassis, we demonstrate that ARCTICA effectively simulates global-scale metabolic flux control. Key findings are that (i) the photosynthetic bioproduction is mainly governed by enzymes within the Calvin–Benson–Bassham (CBB) cycle, rather than by those involve in the biosynthesis of the end-product, (ii) the catalytic capacity of the CBB cycle limits the photosynthetic activity and downstream pathways and (iii) ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) is a major, but not the most, limiting step within the CBB cycle. Predicted metabolic reactions qualitatively align with prior experimental observations, validating our modelling approach. ARCTICA serves as a valuable pipeline for understanding cellular physiology and predicting rate-limiting steps in genome-scale metabolic networks, and thus provides guidance for bioengineering of cyanobacteria.

Genome-scale metabolic network model and phenome of solvent-tolerant Pseudomonas putida S12

BackgroundPseudomonas putida S12 is a gram-negative bacterium renowned for its high tolerance to organic solvents and metabolic versatility, making it attractive for various applications, including bioremediation and the production of aromatic compounds, bioplastics, biofuels, and value-added compounds. However, a metabolic model of S12 has yet to be developed.ResultsIn this study, we present a comprehensive and highly curated genome-scale metabolic network model of S12 (iSH1474), containing 1,474 genes, 1,436 unique metabolites, and 2,938 metabolic reactions. The model was constructed by leveraging existing metabolic models and conducting comparative analyses of genomes and phenomes. Approximately 2,000 different phenotypes were measured for S12 and its closely related KT2440 strain under various nutritional and environmental conditions. These phenotypic data, combined with the reported experimental data, were used to refine and validate the reconstruction. Model predictions quantitatively agreed well with in vivo flux measurements and the batch cultivation of S12, which demonstrated that iSH1474 accurately represents the metabolic capabilities of S12. Furthermore, the model was simulated to investigate the maximum theoretical metabolic capacity of S12 growing on toxic organic solvents.ConclusionsiSH1474 represents a significant advancement in our understanding of the cellular metabolism of P. putida S12. The combined results of metabolic simulation and comparative genome and phenome analyses identified the genetic and metabolic determinants of the characteristic phenotypes of S12. This study could accelerate the development of this versatile organism as an efficient cell factory for various biotechnological applications.

Development and applications of genome-scale metabolic network models.

The genome-scale metabolic network model is an effective tool for characterizing the gene-protein-response relationship in the entire metabolic pathway of an organism. By combining various algorithms, the genome-scale metabolic network model can effectively simulate the influence of a specific environment on the physiological state of cells, optimize the culture conditions of strains, and predict the targets of genetic modification to achieve targeted modification of strains. In this review, we summarize the whole process of model building, sort out the various tools that may be involved in the model building process, and explain the role of various algorithms in model analysis. In addition, we also summarized the application of GSMM in network characteristics, cell phenotypes, metabolic engineering, etc. Finally, we discuss the current challenges facing GSMM.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing.

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Co-Cultivated Enzyme Constraint Metabolic Network Model for Rational Guidance in Constructing Synthetic Consortia to Achieve Optimal Pathway Allocation Prediction.

Synthetic consortia have emerged as a promising biosynthetic platform that offers new opportunities for biosynthesis. Genome-scale metabolic network models (GEMs) with complex constraints are extensively utilized to guide the synthesis in monocultures. However, few methods are currently available to guide the rational construction of synthetic consortia for predicting the optimal allocation strategy of synthetic pathways aimed at enhancing product synthesis. A standardized method to construct the co-cultivated Enzyme Constraint metabolic network model (CulECpy) is proposed, which integrates enzyme constraints and modular interaction scale constraints based on the research concept of "independent + global". This method is applied to construct several synthetic consortia models, which encompassed different target products, strains, synthetic pathways, and compositional structures. Analyzing the model, the optimal pathway allocation and initial inoculum ratio that enhance the synthesis of target products by synthetic consortia are predicted and verified. When comparing with the constructed co-culture synthesis system, the normalized root mean square error of all optimal theoretical yield simulations is found to be less than or equal to 0.25. The analyses and verifications demonstrate that the method CulECpy can guide the rational construction of synthetic consortia systems to facilitate biochemical synthesis.

Genome sequencing and metabolic network reconstruction of a novel sulfur-oxidizing bacterium Acidithiobacillus Ameehan.

Sulfur-oxidizing bacteria play a crucial role in various processes, including mine bioleaching, biodesulfurization, and treatment of sulfur-containing wastewater. Nevertheless, the pathway involved in sulfur oxidation is highly intricate, making it complete comprehension a formidable and protracted undertaking. The mechanisms of sulfur oxidation within the Acidithiobacillus genus, along with the process of energy production, remain areas that necessitate further research and elucidation. In this study, a novel strain of sulfur-oxidizing bacterium, Acidithiobacillus Ameehan, was isolated. Several physiological characteristics of the strain Ameehan were verified and its complete genome sequence was presented in the study. Besides, the first genome-scale metabolic network model (AMEE_WP1377) was reconstructed for Acidithiobacillus Ameehan to gain a comprehensive understanding of the metabolic capacity of the strain.The characteristics of Acidithiobacillus Ameehan included morphological size and an optimal growth temperature range of 37-45°C, as well as an optimal growth pH range of pH 2.0-8.0. The microbe was found to be capable of growth when sulfur and K2O6S4 were supplied as the energy source and electron donor for CO2 fixation. Conversely, it could not utilize Na2S2O3, FeS2, and FeSO4·7H2O as the energy source or electron donor for CO2 fixation, nor could it grow using glucose or yeast extract as a carbon source. Genome annotation revealed that the strain Ameehan possessed a series of sulfur oxidizing genes that enabled it to oxidize elemental sulfur or various reduced inorganic sulfur compounds (RISCs). In addition, the bacterium also possessed carbon fixing genes involved in the incomplete Calvin-Benson-Bassham (CBB) cycle. However, the bacterium lacked the ability to oxidize iron and fix nitrogen. By implementing a constraint-based flux analysis to predict cellular growth in the presence of 71 carbon sources, 88.7% agreement with experimental Biolog data was observed. Five sulfur oxidation pathways were discovered through model simulations. The optimal sulfur oxidation pathway had the highest ATP production rate of 14.81 mmol/gDW/h, NADH/NADPH production rate of 5.76 mmol/gDW/h, consumed 1.575 mmol/gDW/h of CO2, and 1.5 mmol/gDW/h of sulfur. Our findings provide a comprehensive outlook on the most effective cellular metabolic pathways implicated in sulfur oxidation within Acidithiobacillus Ameehan. It suggests that the OMP (outer-membrane proteins) and SQR enzymes (sulfide: quinone oxidoreductase) have a significant impact on the energy production efficiency of sulfur oxidation, which could have potential biotechnological applications.

Elucidating the impact of in vitro cultivation on Nicotiana tabacum metabolism through combined in silico modeling and multiomics analysis.

The systematical characterization and understanding of the metabolic behaviors are the basis of the efficient plant metabolic engineering and synthetic biology. Genome-scale metabolic networks (GSMNs) are indispensable tools for the comprehensive characterization of overall metabolic profile. Here we first constructed a GSMN of tobacco, which is one of the most widely used plant chassis, and then combined the tobacco GSMN and multiomics analysis to systematically elucidate the impact of in-vitro cultivation on the tobacco metabolic network. In-vitro cultivation is a widely used technique for plant cultivation, not only in the field of basic research but also for the rapid propagation of valuable horticultural and pharmaceutical plants. However, the systemic effects of in-vitro cultivation on overall plant metabolism could easily be overlooked and are still poorly understood. We found that in-vitro tobacco showed slower growth, less biomass and suppressed photosynthesis than soil-grown tobacco. Many changes of metabolites and metabolic pathways between in-vitro and soil-grown tobacco plants were identified, which notably revealed a significant increase of the amino acids content under in-vitro condition. The in silico investigation showed that in-vitro tobacco downregulated photosynthesis and primary carbon metabolism, while significantly upregulated the GS/GOGAT cycle, as well as producing more energy and less NADH/NADPH to acclimate in-vitro growth demands. Altogether, the combination of experimental and in silico analyses offers an unprecedented view of tobacco metabolism, with valuable insights into the impact of in-vitro cultivation, enabling more efficient utilization of in-vitro techniques for plant propagation and metabolic engineering.

MetNetComp: Database for Minimal and Maximal Gene-Deletion Strategies for Growth-Coupled Production of Genome-Scale Metabolic Networks.

Growth-coupled production, in which cell growth forces the production of target metabolites, plays an essential role in the production of substances by microorganisms. The strains are first designed using computational simulation and then validated by biological experiments. In the simulations, gene-deletion strategies are often necessary because many metabolites are not produced in the natural state of the microorganisms. However, such information is not available for many metabolites owing to the requirement of heavy computation, especially when many gene deletions are required for genome-scale models. A database for such information will be helpful. However, developing such a database is not straightforward because heavy computation and the existence of replaceable genes render difficulty in efficient enumeration. In this study, the author developed efficient methods for enumerating minimal and maximal gene-deletion strategies and a web-based database system. MetNetComp provides information on 1) a total of 85,611 gene-deletion strategies excluding apparent duplicate counting for replaceable genes for 1,735 target metabolites, 11 constraint-based models, and 10 species; 2) necessary substrates and products in the process; and 3) reaction rates that can be used for visualization. MetNetComp is helpful for strain design and for new research paradigms using machine learning.

Reconstruction and metabolic profiling of the genome-scale metabolic network model of Pseudomonas stutzeri A1501

Pseudomonas stutzeri A1501 is a non-fluorescent denitrifying bacteria that belongs to the gram-negative bacterial group. As a prominent strain in the fields of agriculture and bioengineering, there is still a lack of comprehensive understanding regarding its metabolic capabilities, specifically in terms of central metabolism and substrate utilization. Therefore, further exploration and extensive studies are required to gain a detailed insight into these aspects. This study reconstructed a genome-scale metabolic network model for P. stutzeri A1501 and conducted extensive curations, including correcting energy generation cycles, respiratory chains, and biomass composition. The final model, iQY1018, was successfully developed, covering more genes and reactions and having higher prediction accuracy compared with the previously published model iPB890. The substrate utilization ability of 71 carbon sources was investigated by BIOLOG experiment and was utilized to validate the model quality. The model prediction accuracy of substrate utilization for P. stutzeri A1501 reached 90 %. The model analysis revealed its new ability in central metabolism and predicted that the strain is a suitable chassis for the production of Acetyl CoA-derived products. This work provides an updated, high-quality model of P. stutzeri A1501for further research and will further enhance our understanding of the metabolic capabilities.

Genome-scale reconstruction of the metabolic network in Streptococcus thermophilus S-3 and assess urea metabolism.

Streptococcus thermophilus is an important strain widely used in dairy fermentation, with distinct urea metabolism characteristics compared to other lactic acid bacteria. The conversion of urea by S. thermophilus has been shown to affect the flavor and acidification characteristics of milk. Additionally, urea metabolism has been found to significantly increase the number of cells and reduce cell damage under acidic pH conditions, resulting in higher activity. However, the physiological role of urea metabolism in S. thermophilus has not been fully evaluated. A deep understanding of this metabolic feature is of great significance for its production and application. Genome-scale metabolic network models (GEMs) are effective tools for investigating the metabolic network of organisms using computational biology methods. Constructing an organism-specific GEM can assist us in comprehending its characteristic metabolism at a systemic level. In the present study, we reconstructed a high-quality GEM of S. thermophilus S-3 (iCH492), which contains 492 genes, 608 metabolites and 642 reactions. Growth phenotyping experiments were employed to validate the model both qualitatively and quantitatively, yielding satisfactory predictive accuracy (95.83%), sensitivity (93.33%) and specificity (100%). Subsequently, a systematic evaluation of urea metabolism in S. thermophilus was performed using iCH492. The results showed that urea metabolism reduces intracellular hydrogen ions and creates membrane potential by producing and transporting ammonium ions. This activation of glycolytic fluxes and ATP synthase produces more ATP for biomass synthesis. The regulation of fluxes of reactions involving NAD(P)H by urea metabolism improves redox balance. Model iCH492 represents the most comprehensive knowledge-base of S. thermophilus to date, serving as a potent tool. The evaluation of urea metabolism led to novel insights regarding the role of urease. © 2023 Society of Chemical Industry.

Immune cell-specific and common molecular signatures in rheumatoid arthritis through molecular network approaches

Rheumatoid arthritis (RA) is an autoimmune disorder and common symptom of RA is chronic synovial inflammation. The pathogenesis of RA is not fully understood. Therefore, we aimed to identify underlying common and distinct molecular signatures and pathways among ten types of tissue and cells obtained from patients with RA. In this study, transcriptomic data including synovial tissues, macrophages, blood, T cells, CD4+T cells, CD8+T cells, natural killer T (NKT), cells natural killer (NK) cells, neutrophils, and monocyte cells were analyzed with an integrative and comparative network biology perspective. Each dataset yielded a list of differentially expressed genes as well as a reconstruction of the tissue-specific protein-protein interaction (PPI) network. Molecular signatures were identified by a statistical test using the hypergeometric probability density function by employing the interactions of transcriptional regulators and PPI. Reporter metabolites of each dataset were determined by using genome-scale metabolic networks. It was defined as the common hub proteins, novel molecular signatures, and metabolites in two or more tissue types while immune cell-specific molecular signatures were identified, too. Importantly, miR-155-5p is found as a common miRNA in all tissues. Moreover, NCOA3, PRKDC and miR-3160 might be novel molecular signatures for RA. Our results establish a novel approach for identifying immune cell-specific molecular signatures of RA and provide insights into the role of common tissue-specific genes, miRNAs, TFs, receptors, and reporter metabolites. Experimental research should be used to validate the corresponding genes, miRNAs, and metabolites.

Local flux coordination and global gene expression regulation in metabolic modeling

Genome-scale metabolic networks (GSMs) are fundamental systems biology representations of a cell’s entire set of stoichiometrically balanced reactions. However, such static GSMs do not incorporate the functional organization of metabolic genes and their dynamic regulation (e.g., operons and regulons). Specifically, there are numerous topologically coupled local reactions through which fluxes are coordinated; the global growth state often dynamically regulates many gene expression of metabolic reactions via global transcription factor regulators. Here, we develop a GSM reconstruction method, Decrem, by integrating locally coupled reactions and global transcriptional regulation of metabolism by cell state. Decrem produces predictions of flux and growth rates, which are highly correlated with those experimentally measured in both wild-type and mutants of three model microorganisms Escherichia coli, Saccharomyces cerevisiae, and Bacillus subtilis under various conditions. More importantly, Decrem can also explain the observed growth rates by capturing the experimentally measured flux changes between wild-types and mutants. Overall, by identifying and incorporating locally organized and regulated functional modules into GSMs, Decrem achieves accurate predictions of phenotypes and has broad applications in bioengineering, synthetic biology, and microbial pathology.

Reconciliation and evolution of Penicillium rubens genome-scale metabolic networks-What about specialised metabolism?

In recent years, genome sequencing of filamentous fungi has revealed a high proportion of specialised metabolites with growing pharmaceutical interest. However, detecting such metabolites through in silico genome analysis does not necessarily guarantee their expression under laboratory conditions. However, one plausible strategy for enabling their production lies in modifying the growth conditions. Devising a comprehensive experimental design testing in different culture environments is time-consuming and expensive. Therefore, using in silico modelling as a preliminary step, such as Genome-Scale Metabolic Network (GSMN), represents a promising approach to predicting and understanding the observed specialised metabolite production in a given organism. To address these questions, we reconstructed a new high-quality GSMN for the Penicillium rubens Wisconsin 54-1255 strain, a commonly used model organism. Our reconstruction, iPrub22, adheres to current convention standards and quality criteria, incorporating updated functional annotations, orthology searches with different GSMN templates, data from previous reconstructions, and manual curation steps targeting primary and specialised metabolites. With a MEMOTE score of 74% and a metabolic coverage of 45%, iPrub22 includes 5,192 unique metabolites interconnected by 5,919 reactions, of which 5,033 are supported by at least one genomic sequence. Of the metabolites present in iPrub22, 13% are categorised as belonging to specialised metabolism. While our high-quality GSMN provides a valuable resource for investigating known phenotypes expressed in P. rubens, our analysis identifies bottlenecks related, in particular, to the definition of what is a specialised metabolite, which requires consensus within the scientific community. It also points out the necessity of accessible, standardised and exhaustive databases of specialised metabolites. These questions must be addressed to fully unlock the potential of natural product production in P. rubens and other filamentous fungi. Our work represents a foundational step towards the objective of rationalising the production of natural products through GSMN modelling.